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For the special case of freely evolving Dirac electrons in 1 + 1 dimensions, Feynman checkerboard 
| paths have previously been used to derive Wigner's arrival-time distribution which includes all 

, arrivals. Here, an attempt is made to use these paths to determine the corresponding distribution 

■ of first-arrival times. Simple analytic expressions are obtained for the relevant components of the 

£\j ' first-arrival propagator. These are used to investigate the relative importance of the first-arrival 

contribution to the Wigner arrival-time distribution and of the contribution arising from interference 
between first and later (i.e. second, third, ... ) arrivals. It is found that a distribution of (intrinsic) 
first-arrival times for a Dirac electron cannot in general be consistently defined using checkerboard 
paths, not even approximately in the nonrelativistic regime. 
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I. INTRODUCTION 

i> ' 

In the past decade there has been considerable interest in deriving and understanding arrival-time distributions for 
quantum particles, using a wide variety of approaches (for recent reviews see jl], ^J, [})). Here we focus on an approach 
based on Feynman paths H. 

Yamada and Takagi [| applied the consistent histories approach, with Feynman paths as particle histories, to 
the problem of deriving an intrinsic arrival-time distribution. They considered the special case of a freely evolving 
nonrelativistic quantum particle in one spatial dimension. In |^| they showed that within their approach one cannot 
classify the histories according to the number of times n that a particle arrives at a given spatial point x = X 
during a specified finite time interval because the amplitude for n arrivals is zero for every finite n. Their qualitative 
explanation was that a typical Feynman path, being nondifferentiable in time, intersects X an infinite number of 
times in the given time interval, leaving zero amplitude for any finite value of n. They further suggested that this 
would also be the case for a quantum particle propagating in three spatial dimensions in the presence of a potential. 

The velocity associated with a typical Feynman path for a nonrelativistic electron is infinite at almost ever y point 
on it. This is not the case for a relativistic electron - the velocity associated with a Feynman path ||, [?[] for a 
nonrelativistic electron is infinite at almost every point on it. Hence, such a path will not intersect X an infinite 
number of times in a given finite time interval and the amplitude for n arrivals need not be zero for every finite 
n. This is the primary motivation for the following attempt to derive the first-arrival-time distribution for a freely 
evolving Dirac electron in 1 + 1 dimensions. In Sec II, a checkerboard path derivation || of Wigner's arrival-time 
distribution || which includes all arrivals is sketched, primarily to introduce the basic notation and concepts used in 
the following sections. In Sec. Ill the first-arrival propagator for this special case is derived. Unfortunately, there is 
interference between the n = 1 and n > 1 contributions to the arrival-time distribution so that the distribution of 
intrinsic first-arrival times cannot be consistently defined using this approach. Calculated results for the n = 1 and 
interference contributions are presented in Sec. IV for two simple cases. Concluding remarks are made in Sec. V. 

II. CHECKERBOARD PATHS AND WIGNER'S ARRIVAL-TIME DISTRIBUTION FOR DIRAC 

ELECTRONS 

Consider the 1 + 1 dimensional free-electron Dirac equation in the form 



ihij^>(x,t) = H^(x,t) = -~ihca z ^-^(x,t) - mc 2 a x ^(x,t) (1) 
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with \&(x, i) = ( v E'+(x, t), ^-(x, t)) T a two-component spinor and 



The velocity operator 
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(2) 



§ = =c* a (3) 



has the orthonormal eigenfunctions 



with eigenvalues +c and — c respectively. Hence, 











x + = 


I 


, X- = 





(4) 



*(*,*) = = *+(M)x+ + *-(M)X- (5) 

and the upper and lower components of the wave function ty(x,t) are, respectively, the amplitudes at (x,t) for the 
right-going (+c) and left-going (— c) velocity eigenstates x+ and The four components of the retarded 2x2 
propagator K(x B , tB', %A, £a) = K{B\A) of the 1 + 1 dimensional free-electron Dirac equation (]!]) are labelled by 
velocity directions (+ or — ) at {xa^a) and at (xb,£.b) and are accordingly defined by 

^ + (x Bl t B )\ = f°° dx f K ++ (x B ,t B ;x Al t A ) K + -{x B ,tB]X A ,tA)\ f^+{x Al t A ) 
ty-(x B ,t B ) J y_ oo A l v i ; s:_ + (x B ,t s ;x y!l ,t J 4) if — (x B ,t B ;x A ,t A ) ) \ *_(au, t A ) 

The left subscript, + or — , on K denotes respectively a right-going or left-going arrival at x B at time t B while the 
right subscript, + or — , denotes respectively a right-going or left-going departure from x A at time t A . 

Subtracting the Hermitean conjugate of (p, multiplied from the right by ^>(x, t), from (|j) multiplied from the left 
by the Hermitean conjugate of ^(x, t) gives the continuity equation 

^(x,t)^(x,t) + ^-^(x,t)ca z ^(x,t) = 0. (6) 
at ox 

It is assumed throughout the paper that the parameters of the initial wave function are such that Dirac's original 
identification of 

p(x,t) = 9*(x,t)9(x,t) = |*+(M)| a + \*-(x,t)\ 2 (7) 

and 

J(x,t) = c¥(x,t)a z ®(x,t) = c[\V+(x,t)\ 2 - |*_(x,<)| 2 ] (8) 

with single-electron probability and probability current densities, respectively, is an adequate approximation. 

In Feynman and Hibb's classic book "Quantum Mechanics and Path Integrals" 0] it is stated that the free-electron 
propagator K(x B ,t B ; x A ,t A ) can be constructed from a model in which a particle going from x A at time t A to x B 
at time t B is constrained to move diagonally in space-time at constant speed c in checker fashion (i.e., forward in 
time with spatial increment ±Ax with Ax = cAt = c(t B — t A )/N for each of N equal time steps At > 0)|2(J. Each 
component of the propagator is obtained as the N — > oo limit of the sum over all iV-step checkerboard paths joining 
{x A , t A ) to (x B , t B ), with the first and last steps appropriately fixed, when the weight associated with a path having R 
(noncompulsary) reversals of direction or corners is taken to be (imc 2 At/H) R . Jacobson and Schulman [JtJ regrouped 
the sum-over-paths into a sum-over-R, i.e. 

K 0a (B; A) = i (mc/2h) lim V <f> Pa (R){imc 2 At/K) R (a = ±; /? = ±) , (9) 

R>0 

where <frp a (R) is the number of 0a paths with R noncompulsary reversals. They also evaluated the four checkerboard 
path integrals, obtaining the following closed-form expressions for the components of the propagator: 

K ++ (B:A) =- C XXT Jl (!t)> (10) 
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2A C Iba V A 

K + _(B;A) = K_ + (B;A) =i-J (^) , (12) 

where A c = Ti/mc is the Compton wavelength of the electron, xba = (xb~xa), tBA = (*b ~ ^a) and Iba = ctba with 
tba = — (vba/c) 2 } 1 ^ 2 the proper time for a particle moving with constant velocity vba = xbaI^ba- Jacobson 

and Schulman also determined the number of reversals, Rq, that gives the maximum contribution to the sum in (|^). 
They obtained Rq = Iba/^c and also showed that the sum is dominated by terms having R within Rii?^ 2 of Rq 
The picture that emerges j?], is one in which the particle always moves with speed c and typically travels a 
distance ~ {Iba/tba)^ > A c between reversals of direction; its motion is Brownian with diffusion constant fi/m only 
on scales much larger than this correlation distance. 

Now consider the problem of deriving an expression for the distribution Tl(T;X) of arrival times T at the spatial 
point x = X for an ensemble of Dirac electrons all prepared in the same initial state ^(x, 0). Following Yamada and 
Takagi, it is assumed that the arrival-time distribution for the fictitious particles of the checkerboard model (with 
tA = 0, xb — X and ts — T), should it be well-defined, can be identified with the desired distribution for actual 
electrons. It should be noted that the expression (f?]) for p(x,t) contains no +/— cross terms arising from interference 
between paths arriving at (x,t) with right-going (+) and those with left-going (— ) velocities. Hence, at least for the 
particles of the checkerboard model, the probability density p(x, t) can be decomposed into two contributions, one 
associated only with right-going arrivals at (x, t) and the other only with left-going arrivals: p(x, t) = p + (x,t) + p_ (x, t) 
with p±(x,t) = \^±(x,t)\ 2 . Now, recall that particles following checkerboard paths move only at speed c and that 
the time between reversals in their directions of motion is « {tBA/T~BA)^c/c for those paths that make the dominant 
contribution to the propagator. Hence, for At much less than this correlation time, nearly all of the particles in the 
spatial interval [X — cAt, X] that are right-going at time T — At should arrive at x = X during the time interval 
[T — At,T], giving the dominant contribution to the number of right-going arrivals at X during that time interval. 
Taking the limit At — > dt = + this leads to the prediction that c p+ (X, T) dt right-going particles arrive at x = X 
during the infinitesimal time interval [T — dt,T]. A similar argument applies for left-going particles in [X, X + cAt] 
at time T — At. Hence, the distribution of arrival times T at the spatial point x — X is given by 

n(T;X) = I1 + (T:X) + I1-(T;X) ; U ± (T;X) = \^ ± (X,T)\ 2 j J dtp(X,t) , (13) 

where T max is the maximum arrival time of practical interest. For the special case of free motion in one spatial 
dimension, the general results presented without derivation or discussion by Wigner simplify to ( |l3| ) B . 

III. FIRST-ARRIVAL PROPAGATOR FOR A DIRAC PARTICLE IN 1+1 DIMENSIONS 

The fundamental constants h and c are set to 1 in the analysis of this section but restored in the final results. 

To begin, suppose that the spatial interval of interest is divided into M pieces each of length Ax — xba/M, 
assuming xb > xa- Suppose further that a particular path of N time-steps consists of P spatial steps of length Ax 
to the right and Q to the left so that N = P + Q and M = P — Q. The resulting space-time grid of path segments 
available to the particle is illustrated in Fig. [j]. Denote by K9^(B; A) the component of the propagator (||) associated 
with first arrivals at X = xb at time T — tg- It is constructed from only those paths that reach xb for the first 
time at time ts- To compute we have to count the number of (5a. paths with R noncompulsary reversals for the 
restricted space-time grid shown in Fig. |l|. For simplicity it is assumed that the initial wave function ^(x,tA = 0) 
is sufficiently well localized to the left of xb = X that K+1_(B; A) and K+}_(B;A) are the only components of the 

first-arrival propagator that need to be considered. For xa < xb the components K y(B;A) and K_-(B; A) of 

course have contributions only from multiple-arrival paths. 

A. Counting the number of first arrival paths with R corners 

Computation of $^ (R) involves counting the number of restricted /3a paths with a given number of corners in a 
lattice. A path in the lattice of Fig. |^ - obtained from Fig. [I] by clockwise rotation through 45° and rescaling by a 
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factor of (2 1 / 2 Ax) _1 - is built up of (1, 0) = — > and (0, 1) = f motions. There are no <— and J, motions because the 
paths of Fig. [TJ move only forward in time. Denote by a /-corner a point on the path that is reached by a — » step 
and is left by a f step and by a r-corner a point that is reached by a f step and left by a — + step. Denote by Ri 
the number of corners of type / and by R r the number of corners of type r in a given path. Any path in the u — v 
lattice with given initial and final points (ua,va) and (ub,Vb) can be completely specified by either the coordinates 
of its / corners or by the coordinates of its r corners. Both specifications are needed in the following derivation of the 
first-arrival propagator. 

First, consider the counting problem without any Pa or first-arrival restrictions on the paths and including com- 
pulsary as well as noncompulsary corners. In this simple case, enumeration of the number of paths in a lattice with a 
given number of /- or r-corners is solved in the following manner |l2j . To count the number of paths with Ri /-corners 
one builds two vectors u and v that contain the integer u and w-coordinates of the R t corners of such a path: 



u= (ui,u 2 , ■ ■ ■ ,u Rl ) 

v = («i,i>2,- • • ,Vr,). (14) 



These coordinates satisify the inequalities 



UA + 1 < U\ < U2 ■ ■ ■ < U Rl < U B 

VA < Vi < v 2 ■ ■ ■ < v Rl < Vb - 1, (15) 

where the sequences are strictly ordered. Denote by $™ r (i?;) the number of paths with Ri /-corners, where the label 
" tuc" is a reminder that the paths are without restriction. This number may be evaluated by observing that there are 
ub — ua integers from which to choose the it-coordinates and vb — va from which to choose the ^-coordinates. The 
required number is then 



Similarly, one obtains 



ayW = ( UB R UA )( VB R VA )- (i7) 



Now consider only those paths which do not cross (touching is allowed for the moment) x = xb before t = t B - 
When considering such restricted paths in the u — v lattice (see Fig. |^) it is convenient to choose the bottom-right 
corner of the accessible region as the origin of the u — v coordinate system so that the region below the diagonal u = v 
is out-of-bounds. 

First consider the case in which the paths are specified by the coordinates of their /-corners. This case is the easier 
of the two because /-corners are diagnostic, i.e. a necessary and sufficient condition for a restricted path is that none 
of its /-corners be on the forbidden u > v side of the diagonal. It is thus required to calculate the number of paths 
with Ri /-corners from (ua,va) to (u b ,vb) with ua — a, va — and 115 = vb = b, say, such that vi > it, for all i. 
This number, $>i(Ri; BA), is the number (nq) minus the number of paths with Ri /-corners such that ui > Vi for at 
least one i. To calculate the latter number, take k to be the largest integer such that Uk > Vk and build from jl5| ) 
the following sequences 

ua + 1 = a + 1 <u\ < u 2 - ■ ■ < u fc _i < Vk+i ■ ■ ■ < v Rl < vb - 1 = b - 1 

va = < vi < v 2 ■ ■ ■ < v k < u k < u k +i ■ ■ ■ < u Rl < u B = b . (18) 

The sequences are ordered (as can be checked, see |l2) ) and there is a one-to-one correspondence between the double 
sequences ( |l5| ) and the double sequences ( |l8| ) . The number of all the double sequences (|l^) is 



b-a-l\ I 6+1 
Ri-l \Ri + l 



(19) 



and therefore &i(Ri \ BA) is given by 



Now consider, the computation of the corresponding number of restricted paths from (a, 0) to (6, b) with R r r- 
corners, i.e. $ r (i? r ; BA). A necessary condition for such a path is, of course, that all of its R r r-corners be on the 
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allowed v > u side of the diagonal. However, this is not a sufficient condition because it is possible for the /-corner 
between two consecutive such r-corners to be on the forbidden side of the diagonal. An additional complication is 
that the first /-corner might precede the first r-corner and/or the last /-corner might follow the last r-corner. A simple 
way to include such /-corners in the analysis is to add the end-points (a, 0) = (uq,Uo) an d (b, b) = (u.Rj.+i, UR r +i) 
to the set of R r r-corners to obtain the set of R r + 2 points {(uq,vq), (ui,V\), (ur t ,vr t ), (uR r +i, UR r +i)}. Now, 
for 1 < i < Rr — 1, (iii+i,Vi) is the (diagnostic) /-corner between the consecutive r-corners (m,Vi) and (it-t+i, fi+i)- 
Hence, $ r (i? r ; P^4) is the number (117]) minus the number of paths with R r r-corners such that U{+i > u, for at least 
one i, with < i < R r to allow for paths for which the first and/or last corner is an /-corner. To calculate the latter 
number, take k to be the smallest integer such that Uk+i > Ufe and construct the sequences 

a = u < ui < u 2 < ■■■ < u k < Vk+i < ... < v Rr < v Rr+ i = b 
1 = vq + 1 < vi < V2 < ••• < Vk < Uk+i < •■■ < UR r < u Rr+ i — 1 = 6 — 1. (21) 



The total number of these double sequences is 



and therefore <I> r (i? r ; BA) is given by 



b-a+1 \ fb-1 

Rr J \ Rr 



(22) 



^-('i'JUj-Ct'jCi 1 )' (23) 

Finally, consider the desired [3a first-arrival paths, which are not allowed to touch x = ig before t = ts- The various 
notations A, (xa, 4a) an d «a)2 1 ^ 2 Ax are reserved for the first point of a checkerboard path and B, (xb, ts) and 
(m_b, wb)2 1 / 2 Acc for the last point. Denote by A a , (xa + aAx, tA + At) and (a a , 0)2 1 / 2 Ax the point on a /3a path at 
time t = tA + At and by Bp, (xb — (3Ax, ts — At), (bp, 6^)2 1 ^ 2 Aa; the point at time t — ts — At. 

It is important to note that those paths which may touch but do not cross the diagonal u = v, extending from (0, 0) 
to (bp, bp), do not touch x = xb before t = ts and hence are first-arrival paths. Hence, the above expressions for 
<f>; (Ri ; BA) and $ r (i? r ; BA) with A, B, a and b replaced by A a , Bp, a a and bp, respectively, can be used to evaluate 

and hence Kj^. It is only necessary, for a given choice of a and j3, to determine a a , bp and the relation between 
R and R[ or R r , keeping in mind that R includes only noncompulsary corners. 

For Pa = ++ it is clear from Fig. | that (b+,b+) = (Q,Q), (u B ,v B ) = (Q + 1,Q), (u A ,v A ) = (Q + 1 - P, 0) 
and (a+,0) = (Q + 2 — P, 0). Then a little thought leads to the conclusion that the number 3>^(i?) is equal to 
Q r (R r ; B + A + ) with the identification R = 2R r — 1. Hence, upon replacing a by a + = Q + 2 — P, b by b + = Q, and 



R r by (R+ l)/2 in (23), it follows that 



*+ +{R) ((/+IV2) ((P+l)/2j ((/+IV2) ((i?+l)/2 

/ p-i W 0-1 \ / p-2 \ / Q 

- {(R + l)/2)[(R-l)/2j {(R-l)/2){(R+l)/2 



(24) 



where the identity y ^ ^J = \'kJ~ \ k ij ^ aS ^ een nse( ^- 

For 13a = +- it is clear from Fig. | that (b+,b+) = (Q-1,Q- 1), (u B ,v B ) = (Q,Q - 1), (a_,0) = (Q - P,Q), 
(ua, va) — (Q — P, —1) and $+l(P) is equal to $z(Pj; J3+^4_) with R = 2Ri. Hence, upon replacing a by a_ = Q — P, 



b by 6+ = Q - 1 and i?/ by R/2 in (g0|), it follows that 



* ! U ' _ • i i?/2 1 ) ( %2 X ) ( (if- 2)/2 ) ( (R + 2)/2 ) ' 25 ' 



with ^-(i) = 1. 



B. Evaluation of the propagators 



The first-arrival-time propagators are expressed as 

K^ + (B;A) = ^- Vto ^+l(R)(imAt) 

c odd R>1 
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Ky_(B;A) = Van ^ $\ l) _(R)(imAt) R 



(26) 



even R>0 



C. The case of K 



CO 



In this case we have from (p6[) and using the approximation 



n m /m\, that becomes exact as n — > oo, 



2A r „ AT 

— lim 

„ N^oc 

odd R>1 

This expression may be transformed |22f| to 



lim V (imAt) R [ — 



/ p(H+l)/2 Q(fl-l)/2 p(A-l)/2 Q(«+l 



)/2 



)/2]![(i?-l)/2]! [(fl-l)/2]![(ii+l)/2]! 



— lim V (imAt) fl (P-Q)- \,L„ — - — — . 

2A, Af^oo } ' v l)/2]!([(i?- l)/2]! 



(27) 



. X B A 
! 



lim > 



odd R> ] 



BA 



i V2A c ; [(i?+l)/2]![(i?-l)/2]! 



A J 



l _BA Y^(\k ( l _BA 

2A ^ l J ' 



fc=0 



2 k 



2X C J fc!(fc + l)! 



(28) 



In the last line R has been replaced by 2k + 1 and the limit N —> oo taken. Comparison with the power series 
representation of the Bessel function, 



^) = (^) n E 



k=0 



k\{k + n)V2 



G) 2k , 



immediately gives 



K { A(B;A) 



xba t i Ira 



X c Iba V A 



(29) 



(30) 



D. The case of K<£l 



In this case we start with the expression 

A) — lim Yl ( lwA< ) 

Za c N^oo '■ — ' 

even R >0 

Steps analogous to those above lead to 



R ( P R/2 Q R/2 P( R - 2 )/ 2 Q(«+2)/2 



(R/2)\ {R/2)\ [(R - 2)/2]! [(R + 2)/2}l 



Jn 



''BA 



C - VBA j ( l_BA_ 
C + VBA V A c 



An interesting equality emerges from the above results, namely 



K ++ {B-A) - K__(B;A)=K ( + 1 i(B;A) 



(31) 



(32) 



(33) 



In fact all paths contributing to the K--(B; A) component of the propagator touch the line x — xb at least twice 

and therefore this set of paths is complementary to the one contributing to K++(B; A) in the limit of N — » oo. 
Finally, 



K^\B;A)=K ++ (B;A)-K^(B-A)=K^(B-A) = g^Z^d j^ 

K^-\B;A)=K + ^B;A)-K^(B;A) = - * J 2 

2A C c + v B A \ A c 



(34) 
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IV. INTERFERENCE BETWEEN FIRST AND LATER ARRIVALS 



The decomposition Kp a (B;A) = K^(B;A) 
axrivals of a particle at xb at leads immediately to the corresponding decomposition t^) 



(2 3) 

KnZ '" ' (B;A) according to first and later (second, third, etc.) 



(2,3, 



[x b , is) for the (3 = ± components of the wave function is). Substitution of the latter expression, with 



t A = 0, xb = X and is = T, into the result ( |13| ) for the arrival-time distribution gives 

U(T;X) = U {1] (T;X) + n( 2 ' 3 '-)(T; X) + II^ 1 x 2 > 3 >~)(T; X) . 
For the special case considered here in which the initial probability density p(xa,0) is negligible for xa > X, 



n (1) (T; A) = C 



/ dx A [K ( +i(X, T; x A , 0)*+(x A , 0) + K^(X, T: x A , 0)*-(x A , 0)] 

J CO 



(35) 



(36) 



is the contribution of first arrivals 

i-X 



n (2,3,...)( T;X ) = Q 

+ c 

is the contribution of later-arrivals and 



f dx A [K^'-\X, T; 0)4< + (xa, 0) + K^'-\x, T; x A , 0)*-{x A) 0)] 

J — oo 

g^a [^i 2 ; 3 '- ^, T; 0)*+(x A , 0) + K ( ^--\X, T; xa, 0)4<_(xa, 0)] 



(37) 



n (ix2,3,...) (T;X) = 2C^[( dx A [K^l(X,T;xA,0)^ + (x A ,0) + K { ^_(X,T;xa,0)^^(xa,0)] 

W— oo 

x / c^a [X^'-^X, T; x A , 0)*+(xa, 0) + K^f'-\X, T; x A , 0)tf _(xa, 0)] 

./ — OO 



(38) 



is the contribution due to interference between first and later arrivals. (C is the normalization factor appearing 
in ([H}])). Of particular interest here is the magnitude of the interference contribution relative to the first-arrival 
contribution in the regime vba « c. 

First, however, briefly consider the regime in which vba is so close to c that the correlation distance A c [l — 
(wba/c) 2 ] -1 / 2 for reversal of direction is sufficiently large that for a typical checkerboard path there is insufficient 
time for more than one arrival at xb- To be more quantitative, assume that the initial amplitude f -(sa^a) of the 
— c velocity eigenstate is completely negligible with respect to the initial amplitude ^F^xa^a) of the +c eigenstate 

so that one need consider only K ++ (B;A) and K V {B]A). Also, for (xb^b) = (X,T) and (xa,^a) — (xa,0) 

assume that xba = X — xa is very close to cIba — cT for those values of T for which II(T; A") is nonnegligible 
and for those values of xa for which p(xa,0) is nonnegligible. In this regime, K^ + (B] A) / K ++ {B\ A) w 1 — 5/2 



(2,3, 
++ 



(B; A)/K ++ (B; A) w 6/2 where v B a = (1 - S)c with 6 « 1. In addition, \K-+(B;A)/K ++ (B;A)\ 

Nl/2 



and K 

{Iba/^<^Ba)Jq{IbaJ^c)/ Jx{Iba/^c) with I ba ~ (2S) 1 ^ 2 ctBA- If 8 is sufficiently small that Iba « A c then, using the 

leading term in rt29 ) for n = and for n = 1, |-K" ^(B; A)/ K ++ (B; A)\ « \ c /xba = A C /(X — xa) which is typically 

very much less than 1 for an initial wave packet ^(xajO) that is well-localized away from x = X. Hence, at least 
to the extent that the concepts of single-particle probability and probability current densities are still meaningful in 
the regime in which vba is very close to c, the interference term is very small and II(T; X) = U^(T;X) to a good 
approximation. Strictly speaking, however, for the special case under consideration the first-arrival-time distribution 
is well-defined only in the limit vba — > c. 

Now, consider the nonrelativistic regime. With the definitions ^±(x, t) = 4>±(x, t) exp(— imc 2 t/h), Mx,t) = 
4>+(x,t) + 0-(x, t) and A0(x,i) = 4> + (x,t) — 0-(x, t), the 1 + 1 dimensional free-electron Dirac equation (|l|) can be 
written 



ihd(j)(x,t)/dt — 
ihd&4>{x,t)/dt = 



-ihc dA(f>(x, t)/dx 
-ihc d(p(x, t)/dx . 



2mc 2 (j)(x, t) 



(39) 



If ihd(j)/dt is negligible with respect to 2mc 2 (f> (with c fixed at its actual value, not set equal to infinity) then 
4> can be replaced by —(ifr/2mc)dA(f>/dx in the second equation of (|39|) to obtain the Schrodinger equation 
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ihdA<j)(x,t)/dt 



-(h 2 /2m)d 2 A(j)(x,t)/dx 2 . If one further assumes that \A<f>(x,t) 



» 



i(x,t)| 2 and identi- 
fies 2~ 1 / 2 A(f>(x, t) with the Schrodinger wave function ips{x,t) then the expressions (0) and (ra) immediately lead 
to the desired nonrelativistic expressions, ps(x,t) — \ips{x,t)\ 2 and Js{x,t) — (h/m)'^[tpg(x,t)dips{x,t)/dx] re- 
spectively, for the nonrelativistic probability and probability current densities. Consistent with these consider- 
ations is the following simple choice of initial (t = 0) wave function *S?(xa,0) for the nonrelativistic regime: 
^+(x a,0) = —g^o(xA) and ^-(xajO) = (1 — g 2 ) 1 ^ 2 ^>o(xA) where g is a real constant very close to 2 -1 / 2 (see below) 
and 9q(x) — (27r) _1 / 4 (Ax) -1 / 2 exp[— (2Ax) _2 (x — xo) 2 + ikox] is a minimum-uncertainty-product gaussian with initial 
centroid xq, initial variance Ace, mean wave vector fco and variance Ak = 1/2 Ax. In the numerical calculations pre- 
sented below the constant g is chosen so that the characteristic velocity v(xa, 0) = J(xa, 0)/p(xa, 0) = (2g 2 — l)c = v 
(independent of xa) is equal to v = hk Q /m with vq « c. 
Now, from ([l0|), ( pi] ) and ( |30| ) it immediately follows that 



K^l(B;A)/K, 



and from (O) and (B2) it follows that 



./__(£; A) = 2v B a/{c±vba) ~ 2v BA /c 



(40) 



K + _,_ + (B;A) 



Jq(Iba/K) + [(c - vba)/(c + vba)]J2(Iba/^c) 

Jq(Iba/Xc) + MIbaAc) - (2vba/c)J 2 (Iba/Xc) 
Jo(Iba/K) 



(41) 



In the regime vba/c << 1 under consideration, Iba/^c ~ ctsA/^c = {c/vba){xbaI 'A c ) >> Xba/^c >> 1. Using the 
leading two terms in Hankcl's asymptotic expansions [|l3| of Jq{z) and Ji{z) for large argument z, i.e. 



Hz) « — 

7TZ 



Ja(*) 



1/2 r 



cos [ z 

4 



1/2 r 



5tt 



cos z — 



1 / 7T 

— sin z 

8z V 4 

15 . / 5tt 

— sin z - 

8z V 4 



(N»i) 

(|z|»l), 



(42) 



gives 



-(B;A) 



vba 



cos 



Lad 



cos 



Iba 
Ac 



16x BA 



Iba 
A c 



Ac 



Ac 



(43) 



For Iba/^c = (4«. + 3)7r/4 with n an integer, the right-hand-side of (|43J) is 1. However, in a well-designed arrival- 
time experiment for the wave function under discussion one would arrange that (X — xq) >> Ax so that p(xa,0) 
is completely negligible for xa > X and also that Ax >> A c so that there is negligible probability of generating 
particle-antiparticle pairs. Hence, X c /xba would be extremely small over the important range of xba- Hence, the 
set of xba values where the right-hand-side of ([l3|) is not close to 2vba/c is of small measure and can be ignored 
when considering integrals over xa, provided that vba/c is not itself extremely small (a rough estimate requires that 
vba/c » X c /16ttx B a)- 

Taking into account that — ^_(xa, 0) 
are negligible when \X — Xq\ >> Ax then leads directly to the estimates 



^+(xa, 0) for \v(xa, 0)| << c and that the terms involving K^\_ and Kl\ 



n (1) (T;X) w (l/2)(2v/c) 2 U(T;X) = 2(v/c) 2 H(T; X) 



(44) 



n (ix2,3,...)( T;X ) w 2(l/2)(2u/c)n(T;X) = 2(v/c)U(T;X) (45) 

for the gaussian wave function under consideration. It should be noted that vba = (X — xa)/T has been approximated 
by v which is consistent with (X — xq) >> Ax in the absence of significant wave packet spreading. Fig. |^ and Fig. |?] 
show results for n(T;0), nW(T; 0)/2(v/c) 2 and n( lx2 < 3 >-)(T; 0)/2(v/c) obtained by numerical evaluation of @) to 
( |38| ) for gaussian wave packets with Ak/ko = 0.02 and 0.2 respectively. In the former case the above estimates are 
excellent approximations; in the latter case, even though wave packet spreading is more important, the estimates still 
provide good approximations for the very large differences in overall scale between the three quantities. 
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V. CONCLUDING REMARKS 

In summing up, it is interesting to make a qualitative comparison of the results of the Feynman path and Bohm 
trajectory approaches for investigating arrival times for Dirac electrons. 

In Bohmian mechanics |L4], [l^. [l6| an electron is postulated to be an actually existing point-like particle and an 
accompanying wave which guides its motion. For a Dirac electron in the presence of a potential V(r, t), the time- 
evolution of the guiding wave \&(f, t) is described by the 3D Dirac equation and the trajectory of the point-like 
particle is determined by the equation-of-motion dr(t)/dt = [J(f, t)/p(f, t)]\f—fny It is further postulated that, for 
an ensemble of electrons all prepared in the same initial state ^(r, 0), the probability of such a particle having initial 
position r' ' = r(t = 0) is given by p(f^°\ 0). The various properties stated below for the intrinsic arrival times of the 
point-like particles of Bohm's theory follow readily from the fact that, for agiven initial wave function, trajectories 
with different starting points f (°) never intersect or even touch each other |17fl . 

Now, the expression for the intrinsic ID arrival-time distribution obtained with either approach can be cast in the 
same form as in classical mechanics, namely 

U(T;X) = U + (T;X) + U^(T;X) ; U±(T;X) = ±J±(X,T) j """" dt[J+(X,t) — J-(X,t)] , (46) 

where J + and J_, respectively, are the right-going and left-going components of the probability current density J. 
The decomposition J = J + + J_ is not uniquely defined. The decomposition associated with the fictitious particles 
of the checkerboard path approach, which move only at the speed of light c, is J± = ± c |^±| 2 while that associated 
with the (assumed) actual particles of the Bohm trajectory approach, each of which moves at a variable speed that 
cannot exceed c, is J± — J0[± J] where is the unit step- function. The former decomposition leads to an arrival-time 
distribution H(T;X) proportional to the probability density p(X,T) while the latter leads to one proportional to the 
absolute value of the probability current density, i.e \J(X,T)\. Moreover, unless one or other (or both) of the two 
components of ^(X, t) is zero for T\ < t < T2, there are more arrivals - many more if X2 — T\ is much larger than the 
Jacobson-Schulman correlation time - of the fictitious particles at X during that time interval than there are of the 
supposed actual particles of Bohm's theory. 

Given the probability current density J(x, t) and using the non-crossing property of Bohm trajectories it is straight- 
forward to decompose the intrinsic arrival-time distribution into contributions from first arrivals, from second arrivals, 
etc. with no interference terms between different orders of arrival |is{ ]. In marked contrast to this, the decomposition 
based on Feynman checkerboard paths in general contains a nonzero interference term between first and later (i.e., 
second, third, ...) arrivals so that from the calculation one cannot extract a well-defined intrinsic first-arrival-time 
distribution. In the nonrelativistic regime this interference term can be very large compared to the first-arrival term. 
Because of this and the extremely small correlation length for reversal of direction (« A c which is only about 10~ 3 of 
the diameter of an atom!), suppression of the interference term by decoherence [|l9| within a time interval much less 
than A c /c in duration immediately following the instant of first arrival would be very difficult, if not impossible, in a 
practical arrival-time measurement. Unless this can be achieved, assuming that the first-arrival times of the fictitious 
particles of the checkerboard model are directly relevant to the arrival times measured in a time-of-flight experiment 
on actual electrons is not justified. 
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while the dashed-line path with a = — and j3 = + does contribute. The restricted domain for evaluation of the first-arrival 
propagator is bounded on the right by the vertical dotted line; its shape is shown in grey. 
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